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Abstract 

Background: Anopheles barbirostris is a vector of malaria in Sri Lanka. The taxon exists as a species complex in the 
Southeast Asian region. Previous studies using molecular markers suggest that there are more than 4 distinct clades 
within the An. barbirostris complex in Southeast Asia. The present study characterizes Sri Lankan An. barbirostris 
using mtDNA cytochrome oxidase subunit I (COI) and ribosomal RNA internal transcribed spacer 2 (ITS2) gene 
sequences. 

Findings: DNA was extracted from morphologically identified An. barbirostris specimens from Sri Lanka, the COI 
and ITS2 regions amplified and their sequences analysed by comparison with other GenBank entries. Maximum 
likelihood trees suggested that Sri Lankan An. barbirostris constitute a different molecular type most closely related 
to clade I. 

Conclusions: Considering the uncorrected p distances between the clade I and Sri Lankan specimens it is fair to 
assume that the specimens collected from widely separated locations in Sri Lanka with morphology characteristic 
of An. barbirostris s.l. form a new molecular type with close resemblance to An. barbirostris s.s from Indonesia and 
Thailand. 
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Background 

Anopheles barbirostris van der Wulp is a vector of mal- 
aria in Sri Lanka [1], India and Southeast Asia, and a 
vector of Brugian filariasis in Southeast Asia [2]. Anoph- 
eles barbirostris and 12 related species form the taxo- 
nomically complex Barbirostris Group of malaria vectors 
in this region ([3] and references therein). An. barbirostris 
s.L is a species complex with up to five sibling species 
whose individual distributions are not well established 
[4-7]. Members of this complex are typically zoophillic, 
develop in still or slow moving fresh water such as in 
swamps and rice fields, and are relatively tolerant of or- 
ganic pollution ([3] and references therein). Larvae were 
recently also found in brackish water (salinity of 4-15 
parts per thousand) in domestic wells and mangrove 
swamps in the Jaffna district of Northern Sri Lanka [8]. 
Morphological differentiation of anopheline sibling species 



* Correspondence: gayan156@gmail.com; noble@jfn.ac.lk 
department of Zoology, Faculty of Science, University of Jaffna, Jaffna 
40000, Sri Lanka 

Full list of author information is available at the end of the article 



is difficult and has previously led to significant misiden- 
tification with implications for vector control, e.g. mos- 
quitoes long identified morphologically as Anopheles 
subpictus sibling species B in Sri Lanka were recently 
shown to be Anopheles sundaicus s.l. by analysis of ribo- 
somal RNA gene sequences [9]. 

Other members of the Barbirostris Group, e.g. An. 
reidi and An. barbrumbosus, have been identified in Sri 
Lanka [10] but the five other species of the Barbirostris 
Subgroup, viz. An. campestris, An. donaldi, An. franciscoi, 
An. hodgkini and An. pollicaris [11] have not been re- 
ported in Sri Lanka. The present study was designed to 
characterize morphologically identified An. barbirostris s.l. 
in Sri Lanka using DNA sequences of the cytochrome oxi- 
dase subunit I (COI) of mtDNA and internal transcribed 
spacer 2 (ITS2) of ribosomal RNA in relation to previous 
DNA sequence data from An. barbirostris s.l. collected 
elsewhere in Southeast Asia [6]. 



O© 2014 Gajapathy et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BiolVlGCl C6ntT3l Commons Attribution License (http://creativecommons.Org/licenses/by/4.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain 
Dedication waiver (http://creativecommons.Org/publicdomain/zero/1.0/) applies to the data made available in this article, 
unless otherwise stated. 



Gajapathy et al. Parasites & Vectors 2014, 7:348 
http://www.parasitesandvectors.eom/content/7/1/348 



Page 2 of 5 



Methods 

Mosquito collection 

Anopheline mosquitoes were collected from the widely 
separated locations of Hambantota, Nainadeevu and 
Puttalam in Sri Lanka (Figure 1) from January to July, 
2012 using cattle baited traps. Mosquitoes were identi- 
fied based on their morphology using taxonomic keys 
created for Sri Lankan anopheline mosquitoes [12]. 

DNA extraction and PCR amplification of COI and ITS2 
regions 

Among 27 morphologically identified An. barbirostris s. 
/., DNA was extracted from two, one and two specimens 
collected from Hambantota, Nainadeevu and Puttalam 
respectively using previously reported methods [13]. The 
COI region was PCR amplified with forward 5' TTG 
ATT TTT TGG TCA TCC AGA AGT 3' and reverse 5' 
TAG AGC TTA AAT TCA TTG CAC TAA TC 3' 
primers [6]. Each PCR (50 uL total volume) reaction 
mixture contained lx Thermo Scientific Ready Mix 
(1.5 mM Mg Cl 2 , 0.625U Taq DNA polymerase, 0.2 mM 
dNTPs and lx Taq buffer), 2-4 ul of template DNA and 
lOpmoles of each primer. The thermal profile for PCR 
was 94°C for 5 min, followed by 35 cycles of 95°C for 




Figure 1 Map of Sri Lanka showing Anopheles barbirostris 
collection sites. 
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40 s, 50°C fori min and72°C fori min. A final extension 
temperature of 72°C was set for 10 min. 

The ITS2 region was amplified using forward 5' ATC 
ACT CGG CTC ATG GAT CG 3' and reverse 5' ATG 
CTT AAA TTT AGG GGG TAG TC 3' primers [6]. 
Fifty ul reaction volumes were prepared with lx reaction 
buffer, 2.5 mM MgCl 2 ; 0.2 mM dNTPs, 20pmoles of 
each primer, 1.0 U of Taq polymerase and 5 ul of tem- 
plate DNA. The PCR thermal profile was 94°C for 
5 min, followed by 35 cycles of 94°C for 1 min, 55°C for 
2 min, 72°C for 2 min and a final extension at 72°C for 
10 min. The amplified products were separated by agar- 
ose gel electrophoresis, stained with ethidium bromide 
and sequenced at Macrogen, Amsterdam, Netherlands. 

Sequences were edited in Finch TV (Geospiza, Seattle, 
USA) and aligned with Clustal W (MEGA 5.1) [14]. A 
phylogenetic tree of the sequences was created using the 
maximum likelihood approach (Phyml 3.0) [15]. A HKY 
substitution model was selected with a gamma distribu- 
tion based on the lower Bayesian criterion index (MEGA 
5.1) [14]. Genbank deposited sequences from Thailand 
(16 for COI and 28 for ITS2) and Indonesia (08 samples 
for COI and 06 samples for ITS2) were also used in the 
tree construction [5,6]. 

Findings 

Anopheles barbirostris specimens were identified based 
on standard morphological characteristics [12]. The 
presence of a median tuft was observed in all flies. Me- 
dian scales in the abdominal sternum were also ob- 
served. Maxillary palps were dark-scaled and the wings 
had a narrow apical fringe spot at the vein R 4+5 . 

All the chromatograms were clean and the ambiguous 
ends were removed before sequence analysis. The phylo- 
genetic tree (Figure 2) constructed with 400 bp of COI 
sequence data, including those from previously pub- 
lished sequences, confirmed that all five sequenced Sri 
Lankan specimens were members of the Barbirostris 
Complex. Sri Lankan specimens formed a sister clade to 
clade I of Paredes-Esquivel et al. [6], which is An. bar- 
birostris s.s. Compared to sequences from the most 
closely related clade I, the Sri Lankan specimens showed 
5.1% variation (with 20 variable sites) in the COI se- 
quences. All the mutations were silent except two where 
tyrosine was replaced by serine and alanine replaced by 
proline. Intra-specific variations (0.4%), not associated 
with specific location, were observed among the Sri 
Lankan COI sequences. 

The maximum likelihood tree (Figure 3) created with 
580 bp of ITS2 sequence data varied slightly from the 
COI tree but the Sri Lankan specimens remained most 
closely related to clade I. Compared to sequences from 
the most closely related clade I, the Sri Lankan speci- 
mens showed 5.4% variation (with 28 variable sites 
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Figure 2 Maximum Likelihood tree created from C01 sequence data (400 bp) using Phyml [15] with the HKY + G + 1 substitution 
model. The nodes are shown with their bootstrap values; An. coustani was used as the out-group. The Genbank accession numbers for Sri Lankan 
samples are KJ580501 - KJ580505. 
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Figure 3 Maximum likelihood tree created from ITS2 sequence data (580 bp) using Phyml [15] with the HKY + G + 1 substitution model. 

The nodes are shown with their bootstrap values; An. coustani was used as the out-group. The GenBank accession numbers for Sri Lankan 
samples are KJ580496 - KJ580500. 
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including 6 indels). The ITS2 sequences from all five Sri 
Lankan specimens were identical. 

All the nodes in both trees (Figures 2 and 3) were well 
supported with strong bootstrap values. The maximum 
likelihood trees created with HKY + G + 1 model in 
Phyml were in agreement with the Bayesian trees (data 
not shown). 

Anopheles barbirostris were previously classified into 
five distinct clades by Paredes-Esquivel et al. [6], based on 
COI and ITS2 data. Clade IV was identified as a species 
with intermediate morphology between An. barbirostris 
and An. campestris. Later, clade I was identified unequivo- 
cally as An. barbirostris s.s. and Clade II was described 
as a new species named An. vanderwulpi [16]. Saeung 
et al. [5] recognized three sibling species in the An. 
barbirostris complex based on karyotyping and COI, 
COII and ITS2 sequences. The sibling species desig- 
nated Al and A2 by Saeung et al. [5] correspond to 
clades III and IV, respectively, of Paredes-Esquivel et 
al. [6]. Differential susceptibility of the members of the 
Barbirostris Complex to infection with malaria para- 
sites has been established in laboratory reared colonies 
in Thailand [17]. 

The disparity between mitochondrial (COI) and nu- 
clear gene (ITS2) trees in the present analysis may be 
explained in part by different evolutionary rates. Mito- 
chondrial DNA typically mutates at a higher rate than 
nuclear sequences [18] such as the ITS2, and thus the 
latter may be more appropriate for analysing relation- 
ships between species. 

Conclusions 

Considering the uncorrected p distances of 4.5% (COI) 
and 5.1% (ITS2) between the clade I and Sri Lankan 
specimens it is fair to assume that the specimens col- 
lected from widely separated locations in Sri Lanka with 
morphology characteristic of An. barbirostris s.l. form a 
new molecular type with close resemblance to the clade 
I {An. barbirostris s.s) of Paredes-Esquivel et al. [6]. 
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